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Abstract 

We use the scale of neutrino mass to derive model-independent naturalness constraints on possible 
contributions to muon decay Michel parameters from new physics above the electroweak symmetry- 
breaking scale. Focusing on Dirac neutrinos, we obtain a complete basis of effective dimension four 
and dimension six operators that are invariant under the gauge symmetry of the Standard Model 
and that contribute to both muon decay and neutrino mass. We show that - in the absence of fine 
F ^ . tuning - the most stringent bounds on chirality-changing operators relevant to muon decay arise 

Ph ! from one-loop contributions to neutrino mass. The bounds we obtain on their contributions to the 

Michel parameters are four or more orders of magnitude stronger than bounds previously obtained 
in the literature. We also show that there exist chirality-changing operators that contribute to 
muon decay but whose flavor structure allows them to evade neutrino mass naturalness bounds. 
We discuss the implications of our analysis for the interpretation of muon decay experiments. 
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I. INTRODUCTION 



Precision studies of muon decay continue to play an important role in testing the Standard 
Model (SM) and searching for physics beyond it. In the gauge sector of the SM, the Fermi 
constant G M that characterizes the strength of the low-energy, four-lepton /x-decay operator 
is determined from the p lifetime and gives one of the three most precisely-known inputs 
into the theory. Analyses of the spectral shape, angular distribution, and polarization of 
the decay electrons (or positrons) probe for contributions from operators that deviate from 
the (V — A) <g) (V — A) structure of the SM decay operator. In the absence of time-reversal 
(T) violating interactions, there exist seven independent parameters - the so-called Michel 
parameters 0, 0| - that characterize the final state charged leptons: two (p, 77) that describe 
the spatially isotropic component of the lepton spectrum; two (£, 5) that characterize the 
spatially anisotropic distribution; and three additional quantities (£', 77") that are needed 
to describe the lepton's transverse and longitudinal polarization 1 . Two additional parame- 
ters {a'/A, (3' /A) characterize a T-odd correlation between the final state lepton spin and 
momenta with the muon polarization: S e ■ k e x S^. 

Recently, new experimental efforts have been devoted to more precise determinations 
of these parameters. The TWIST Collaboration has measured p and 5 at TRIUMF0, 0], 
improving the uncertainty over previously reported values by factors of ~ 2.5 and ~ 3, 
respectively. An experiment to measure the transverse positron polarization has been carried 
out at the Paul Scherrer Institute (PSI), leading to similar improvements in sensitivity 
over the results of earlier measurements [5J]. A new determination of P M £ with a similar 
degree of improved precision is expected from the TWIST Collaboration, and one anticipates 
additional reductions in the uncertainties in p and 5\^. 

At present, there exists no evidence for deviations from SM predictions for the Michel 
parameters (MPs). It is interesting, nevertheless, to ask what constraints these new mea- 
surements can provide on possible contributions from physics beyond the SM. It has been 
conventional to characterize these contributions in terms of a set of ten four-fermion opera- 



where the sum runs over Dirac matrices T 7 = 1 (S), 7° (V), and cr Q/3 /v^2 (T) and the 
subscripts p and e denote the chirality (R,L) of the muon and final state lepton, respectively 2 . 
In the SM, one has g\r = 1 and all other = 0. A recent, global analysis by Gagliardi, 
Tribble, and Williams [q] give the present experimental bounds on the g] that include the 
impact of the latest TRIUMF and PSI measurements. 

Theoretically, the g] can be generated in different scenarios for physics beyond the SM. 
The most commonly cited illustration is the minimal left-right symmetric model that gives 
rise to non-zero g RR , g RL , and g\ R . From a model-independent standpoint, the authors of 
Ref. j^l recently observed that the operators in Eq. having different chiralities for the 
muon and final state charged lepton will also contribute to the neutrino mass matrix m^ B 
through radiative corrections. Consequently, one expects that the present upper bounds on 



1 The parameters r\ and rf' are alternately written in terms of the independent parameters a/ A and (3 /A. 

2 The normalization of the tensor terms corresponds to the convention adopted in Ref. |jj . We do not specify 
the neutrino flavors in Eq. since the /u-decay experiments do not observe the final state neutrinos. 
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m v should imply bounds on the magnitudes of the gl . The authors of Ref. [9| argued that the 
most stringent limits arise from two-loop contributions because the one-loop contributions 
are suppressed by three powers of the tiny, charged lepton Yukawa couplings. The two- 
loop constraints are nonetheless stronger than the present bounds give in Ref. and could 
become even more so with the advent of future terrestrial and cosmological probes of the 
neutrino mass scale. 

In this paper, we present the results of a follow-up analysis of m v constraints on the 
//-decay parameters, motivated by the observations of Ref. [9| and the new experimental 
developments in the field. Our study follows the approach of Ref. (HI U| used recently in 
deriving model-independent naturalness bounds on neutrino magnetic moments implied by 
the scale of m v . For concreteness, we work with an effective theory that is valid below a 
scale A lying above the weak scale v ~ 246 GeV and that contains SU(2)ixU(l)y-invariant 
operators built from Standard Model fields plus right-handed (RH) Dirac neutrinos 3 . We 
consider all relevant operators up to dimension n = 6 that could be generated by physics 
above the scale A. For simplicity, we restrict our attention to two generations of lepton 
doublets and RH neutrinos. Extending the analysis to include a third generation increases 
the number of relevant operators but does not change the substantive conclusions. While the 
spirit of our work is similar to that of Ref. the specifics of our analysis and conclusions 
differ in several respects: 

i) The effective theory that we adopt allows us to compute contributions to m u from 
scales lying between the weak scale v and the scale of new physics A. In contrast, 
the authors of Ref. used a Fierz transformed version of £M" deca y in Eq. (JU), which 
is not invariant under the SM gauge group and, therefore, should be used to analyze 
only contributions below the weak scale. 

ii) We show that for the two flavor case the operators in £^- deca y proportional to g^ R 
and g RL arise from twelve independent dimension n = 6 gauge-invariant four-fermion 
operators, while those containing g\ R and g RL are generated by four independent n = 6 
operators that contain two fermions and two Higgs scalars. 

iii) While the operators that contribute to /x-decay have dimension n = 6 or higher, the 
lowest dimension neutrino mass operator occurs at n = 4. The authors of Ref. used 
dimensional regularization (DR) to analyze the mixing between the n = 6 yU-decay 
and neutrino mass operators, but did not consider mixing with the n = 4 operator 
that cannot be determined with DR. We derive order-of-magnitude constraints on the 
n = 6 operator coefficients implied by this mixing, which depends only linearly on the 
lepton Yukawa couplings and which gives the dominant constraints for A ^> v. 

iv) For A not too different from v, constraints associated with mixing among the n = 6 
operators can, in principle, be comparable to those arising from contributions to the 
n = 4 mass operator. We carry out a complete, one-loop analysis of this mixing and 
show that only the neutrino magnetic moment and two-fermion/two-Higgs operators 
mix with the n = 6 neutrino mass operator to linear order in the lepton Yukawa 
couplings. We derive the resulting bounds on the g\ R RL that follow from this mixing 



We defer a study of the effective theory containing Majorana neutrinos to a subsequent publication. 
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and find that they are comparable to those associated with induced n = 4 mass 
operator for A > v. 

v) From the mixing with the n = 4 and n = 6 mass operators, we find that the bounds 
on the are two or more orders of magnitude stronger than obtained in Ref. J9| 
and at least three orders of magnitude below the experimental limits given in Ref. |S|] . 

vi) The neutrino mass implications for the couplings g^ RRL are more subtle. Of the twelve 
independent four-fermion operators that contribute to these couplings, only eight are 
directly constrained by the scale of neutrino mass and naturalness considerations. 
Their contributions to the Q^rrl are generally ~ 10 4 times smaller than the present 
experimental bounds, and ~ 10 3 times smaller than obtained in the analysis of Ref. 0. 
We show, however, that the flavor structure of the remaining four operators allows 
them to evade such constraints. While from a theoretical perspective one might not 
expect their contributions to be substantially larger than those from the constrained 

S T 

operators, experimental efforts to determine the g LR RL remain a worthwhile endeavor. 

A summary of our results is given in Table |1J In the remainder of the paper we give 
the details of our analysis. In Section m we write down the complete set of independent 
operators through n = 6 that contribute to m AB and/or /i-decay. Section ITTT1 gives our 
analysis of operator mixing, while in Section IIVI we discuss the resulting constraints on the 
9lrrl that follow from this analysis and the present upper bounds on the neutrino mass 
scale. We summarize in Section |3 The Appendix discusses the sensitivity of the Michel 
parameters to the neutrino magnetic moment operators. 



II. OPERATOR BASIS 

To set notation, we follow Ref. jib] and consider the effective Lagrangian 

where fi is the renormalization scale, n > 4 is the operator dimension, and j is an index 
running over all independent operators of a given dimension. The lowest dimension neutrino 
mass operator is 

0$ AD = l A ^E (3) 

where L A is the left-handed (LH) lepton doublet for generation A, u R is a RH neutrino 
for generation D and </> = iT2<f)* with <ft being the Higgs doublet field. After spontaneous 
symmetry breaking, one has 



so that 







r 4 ro (4) -> -m AD v A v D 

AD 



(4) 



m AU = ~C* AD v/V2 . (5) 
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TABLE I: Constraints on /i-decay couplings gja- The first eight rows give naturalness bounds in 
units of (v/A) 2 x (m^/leV) on contributions from n = 6 muon decay operators (defined in Section 
[U] below) based on one-loop mixing with the n = 4 neutrino mass operators. The ninth row gives 
upper bounds derived from a recent global analysis of Ref. 0] , while the last row gives estimated 
bounds from Ref. derived from two-loop mixing of n = 6 muon decay and mass operators. A 
"-" indicates that the operator does not contribute to the given g]^, while "None" indicates that 
the operator gives a contribution unconstrained by neutrino mass. The subscript D runs over the 
two generations of RH Dirac neutrinos. 

Source \gfj \gj R \ bjj j^j jg^ \g^\ 

4 x i(r 7 2 x i(r 7 - 

4 x i(r 7 ----- 

None None - 

8 x 1(T 5 4 x KT 5 
8 x KT 5 
None None 

8 x KT 7 

2 x KT 4 



Global [g] 0.088 0.025 0.417 0.104 0.036 0.104 
Two-loop [9] 10" 4 10" 4 10~ 2 10~ 2 10" 4 10~ 2 

The other n = 4 operators are those of the SM and we do not write them down explicitly 
here. 

For the case of Dirac neutrinos that we consider here, there exist no gauge-invariant n = 5 
operators. In considering those with dimension six, it is useful to group them according to 
the number of fermion, Higgs, and gauge boson fields that enter: 

Four fermion: 



q(6) 

F, 122D 
(6) 

F, 212D 
(6) 

F, 11 2D 

W F,211D 
(6) 

F, 121D 
(6) 

F, 221D 
G (6) 

V,2D 
(6) 

V,1D 



Vr^VrVr^^Vr 

L£ R I R L 

LvrVrL 

e ij Li£ R ZjV R 

Here £r is the right-handed charged lepton field. Several of the operators appearing in this 
list can contribute to /x-decay, but only the last one can also contribute to m^ D through 
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radiative corrections. Including flavor indices, we refer to this operator as 

F, ABCD — e l R- L j V R K°) 

where the indices i,j refer to the weak isospin components of the LH doublet fields and 
Fermion-Higgs: 



t{L A Yr a L B ){<P + T a D,<P) 

i{ihn B R ){^D,cj>) (7) 
i{9irv B ){<t> + D,<j>) 

Neither of the first two operators in the list (J7J) can contribute significantly to m AD since they 
contain no RH neutrino fields. Any loop graph through which they radiatively induce m AD 
would have to contain operators that contain both LH and RH fields, such as 0$ AB or other 
n = 6 operators. In either case, the resulting constraints on the operator coefficients will be 
weak. For similar reasons, the third and fourth operators cannot contribute substantially 
because they contain an even number of neutrino fields having the same chirality and since 
the neutrino mass operator contains one LH and one RH neutrino field. Only the last 
operator 

Of^iili^i^D^) (8) 

can contribute signficantly to m„ since it contains a single RH neutrino. It also contributes 
to the //-decay amplitude after SSB via the graph of Fig. since the covariant derivative 
contains charged W-boson fields. We also write down the n = 6 neutrino mass operators 

0^ AD = {L A 4>v°){<P+<P) (9) 

as well as the charged lepton mass operator (L<f)£ii)(<j) + <j)) that we do not use in the present 
analysis. 

Fermion-Higgs- Gauge: 
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Lr a ^D u LW a uv 



L^D U LB. 



UrYD^VrB^ (10) 
gi {La" v (j))lRB„ v 

9l (La^)u R B^ 

As for the fermion-Higgs operators, the operators in (JIUj) that contain an even number of vr 
fields will not contribute significantly to m AB , so only the last two in the list are relevant: 

Of, ad = gi(L A a^)v°B, u (11) 
<D$ AD = g 2 (L A a^ry)^W^ (12) 

In addition to these operators, there exist additional n = 6 operators that contain two 
derivatives. However, as discussed in Ref. [l(|, they can either be related to Of AD and 

Ow ad through the equations of motion or contain derivatives acting on the vr fields so 
that they do not contribute to the neutrino mass operator. Consequently, we need not 
consider them here. We also observe that the operator OyJ AD will also contribute to the 
/x-decay amplitude via graphs as in Fig. [TJd. We have computed its contributions to the 
Michel parameters and find that they are suppressed by (~ m M /A) 2 < 1.7 x 10 -7 relative 
to the effects of the other n = 6 operators. This suppression arises from the presence of 
the derivative acting on the gauge field and the absence of an interference between the 
corresponding amplitude and that of the SM. Finally, we note that the operators whose 
chiral structure suppresses their contributions to the neutrino mass operator (as discussed 
above) may, in general, contribute to muon decay via the terms in Eq. having e = fi. We 
do not consider these terms in this study. 




FIG. 1: Contributions from the operators (a) O-n- * D an d 0>) ®w AD (denoted by the shaded box) 
to the amplitude for fj,- decay. Solid, dashed, and wavy lines denote fermions, Higgs scalars, and 
gauge bosons, respectively. After SSB, the neutral Higgs field is replaced by its vev, yielding a 
four-fermion /i-decay amplitude. 
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III. OPERATOR MIXING 



In analyzing mixing among operators that contribute to both /z-decay and m A it is 
useful to consider separately two cases: (i) mixing between the n = 6 operators that enter 
/x-decay and the n = 4 mass operator, AD , and (ii) mixing among the relevant n = 6 
operators. In general, contributions to m AD involving the second case will be smaller than 
those that involve mixing with AD by ~ (v/A) 2 , since Ofj AD contains an additional 
factor of (<f>* '(f>) I A 2 . We first consider this case and employ dimensional analysis to derive 
neutrino mass naturalness bounds on the n = 6 operator coefficients. For v not too different 
from A, the impact of the n = 6 mixing can also be important, and in this case we can 
employ a full renormalization group (RG) analysis to derive robust naturalness bounds. 



A. Mixing with 0$ AD 

The analysis of Ref. employed dimensional regularization (DR) to regularize the one- 
and two-loop graphs through which four-fermion operators containing a single z/r field con- 
tribute to the n = 6 mass operator. Mixing with lower-dimension operators cannot be 
treated using DR since the relevant graphs are quadratically divergent and must be pro- 
portional to the square of a mass scale. For fi > v, all fields are massless, and /i itself 
appears only logarithmically. Since the mass operator exists for zero external momentum, 
all quadratically-divergent graphs vanish in this case. 

The n = 4 mass operator will nevertheless receive contributions from the quadratically 
divergent graphs containing the n = 6 operators. Since the integration region includes 
momenta of order the cut-off, these quadratically divergent graphs will have magnitude 
~ A 2 / (47r) 2 . In DR, these contributions are absorbed into the operator coefficients AD 
whose values are taken from experiment. One may, however, estimate the size of these 
contributions either using a gauge-invariant regulator, such as the generalized Pauli-Villars 
regulator of Ref. 0], or using dimensional arguments. Since we are interested in order-of- 
magnitude constraints, use of the latter is sufficient. 

The relevant one-loop graphs are shown in Fig. |21 For the mixing of the four-fermion 
operators Op ABCD into 0$ AD , two topologies are possible, associated with either the fields 
(L A , is %) or (L B , Vft) living on the external lines. For the mixing of Op\ BCD as well as of 
C| AB into 0^ AEn one insertion of the Yukawa interaction f Ac l B L A is needed to convert 
the internal, RH lepton into a LH one. In contrast, no Yukawa insertion is required for the 
mixing of { B ] AD and 0$ AD into Of l AD . 

To simplify the analysis of mixing with the Op\ BCD we note that one may always redefine 
the fields L A and £ B so that the charged lepton Yukawa matrix Jad is diagonal. Specifically, 
we take 

L A - L Al = S AB L B (13) 
£ c ^ f< = T CD l D 



with Sab and Tqd chosen so that 



Lf£ = L'f diag £' (14) 
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FIG. 2: One loop graphs for the mixing of the n = 6 operators (denoted by the shaded box) into the 
n = 4 mass operator AD . Solid, dashed, and wavy lines denote fermions, Higgs scalars, and 
gauge bosons, respectively. Panels (a,b,c) illustrate mixing of 0^ w , and Op, respectively, 

into { m] ad - 

where L, L' denote vectors in flavor space, / denotes the Yukawa matrix in the original 
basis, and /diag = f T. We note that the field redefinition (|T3j) differs from the conven- 
tional flavor rotation used for quarks, since we have performed identical rotations on both 
isospin components of the left-handed doublet. Consequently, gauge interactions in the new 
basis entail no transitions between generations. We also note that Eqs. (|13|) also imply a 
redefinition of the operator coefficients C\ f AD , C F ABCD , etc.. For example, one has 

U M,A'D — M, AD "AT, A' A l-LOJ 

s~ifii /^y6 c C T* 

^F.A'B'C'D — ^F,ABCD °A'A>~>B'B J-C'O 



where a sum over repeated indices is implied. Diagonalization of the neutrino mass matrix 
requires additional, independent rotations of the fields after inclusion of radiative con- 
tributions to the coefficients C^f AD generated by physics above the weak scale. Since we 
are concerned only with contributions generated above the scale of SSB, we will not perform 
the latter diagonalization and carry out computations using the L', £' R basis 4 . 

In this case, the only four fermion operators F ABCD that can contribute substantially 
to m AD are those having either A = C or B = C. Thus, we obtain the following estimates 



For notational simplicity, we henceforth omit the prime superscripts. 
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of the contributions from the n = 6 operators to the coefficient of the n = 4 mass operator: 



O 



(6) 



B, AD 



o 



(6) 

W,AD 



O 



0(6) 
V,AD 

(6) 

F, ABAD 







(6) 

F, ABBD 



/-»4 

°M, AD 



U M AD 



^4 

°M, AD 



^4 

°Jlf, BD 



^4 

U M, AD 



c 



47T COS 2 9w 

3a 
47r sin 2 

JAA „6 
16n 2 U V,AD 

f. 



B, AD 



/~f6 

^W,AD 



(16) 



AA 



4tt 2 
/bb 



^F, ABAD 



167T 



C6 
- F, ABBD 



,(6) 



where (9y^ is the weak mixing angle. 

The relative factor of 3 cot 2 6 W for the mixing of 0^ AD compared to the mixing of <D% AD 
arises from the ratio of gauge couplings {g/g') 2 and the presence of a r • f appearing in Fig. 

The factor of two that enters the mixing of Op\ BAD compared to that of O f \ bbd 
arises from the trace associated with the closed chiral fermion loop that does not arise for 

^F,ABBD- 

We observe that there exist two four-fermion operators that contribute to /z-decay that 
do not contribute to C\ t AD in the basis giving a diagonal Jab'- ®faabd with either A = 
1, B = 2 or A = 2, B = 1. As we discuss in Section HVl these operators contribute to g£^ 

S T 

and g BL , respectively. Consequently, the magnitudes of these couplings are not directly 
bounded by m u and naturalness considerations, as indicated in Table HJ From a theoretical 
standpoint, one might expect the magnitudes of C F ll2D and C F22XD to be comparable to 
those of the other four-fermion operator coefficients in models that are consistent with the 
scale of neutrino mass. Nevertheless, we cannot a priori rule out order of magnitude or more 
differences between operator coefficients. 



B. Mixing among n = 6 operators 

Because Of I AD contains one power of (<f)^<f))/A 2 compared to 0$ AD , the constraints 
obtained from mixing with the former will generally be weaker by ~ (v/A) 2 . However, for 
A not too different from the weak scale, the n = 6 mixing can be of comparable importance 
to the n = 4 case. Here, we study the mixing among n = 6 operators by computing all 
one-loop graphs that contribute using DR and performing a renormalization group (RG) 
analysis. Doing so provides the exact result for contributions to the one-loop mixing from 
scales between A and v, summed to all orders in f AA ln(v/A) and aln(u/A). 

In carrying out this analysis, it is necessary to identify a basis of operators that close under 
renormalization. We find that the minimal set consists of seven operators that contribute 
to /i-decay and m^ D : 

ro( 6 ) h7\ 

W B,AD' ^W,ADi W M,AD' W v,AD> F AAAD > F, ABBD ' W F,BABD • 

For simplicity, we have included a single RH neutrino field u B in all seven operators. While 
one could, in principle, allow for different vr generation indices, the essential physics can be 
extracted from an analysis of this minimal basis. 
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The classes of graphs relevant to mixing among these operators are illustrated in Fig. 
El where we show representative contributions to operator self-renormalization and mixing 
among the various operators. The latter include mixing of all operators into Ofj AD (a-c); 

mixing of 0^ AD , 0^ AD , and 0^ AD into 0^ AD (d,e); and mixing between the four-fermion 

operators and the magnetic moment operators (f,gL Representative self-renormalization 
graphs are given in Fig. [HJ^h-j) . As noted in Ref. |9|, the mixing of the the four-fermion 

operators into Ofl AD contains three powers of the lepton Yukawa couplings and is highly 
suppressed. In contrast, all other mixing contains at most one Yukawa insertion. 

Working to first order in the Jaa we find a total of 59 graphs that must be computed, 
not including wavefunction renormalization gra phs that are not shown. Twenty-two of these 



a pns tf. 

graphs were computed by the authors of Ref. [1QJ in their analysis of the mixing between 
O^r a n and the magnetic moment operators. Here, we compute the remaining 37. As in 
Ref. [10], we work with the background field gauge [l3j] in d = 4 — 2e spacetime dimensions. 
We renormalize the operators using minimal subtraction, wherein counterterms simply re- 
move the divergent, 1/e terms from the one-loop amplitudes. The resulting renormalized 
operators Of^ are expressed in terms of the unrenormalized operators O^f 1 as 

= E ZjiZF^O? = £ ZjiO® , (18) 

k k 

where 

(6) = z n L /2 z n,,2 ^) (lg) 

are the //-independent bare operators, Z^ 2 and Z^j 2 are the wavefunction renormalization 
constants for the fields L A and 0, respectively; ni and are the number of LH lepton and 

Higgs fields appearing in a given operator; and Z~^Z r ^^ 2 Z n ^ 2 are the counterterms that 
remove the 1/e divergences. 

Since the bare operators do not depend on the renormalization scale, whereas the 

Zjj} and the Of^ do, the operator coefficients C® must carry a compensating //-dependence 
to ensure that £ c fr is independent of scale. This requirement leads to the RG equation for 
the operator coefficients: 

+ E °l 7* = ( 2 °) 

where 
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is the anomalous dimension matrix. We obtain 5 
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27q 2 /^ 



3ai 

47T 



9ai/AA 

47T 
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3(3ai— 02) 
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9ai/.B 
2tt 

9Q2/-B 
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9a 1/. 



BB 



4tt 
9a2/Bfl 
4tt 



3(ai+a 2 ) 3(ai-a 2 ) 
8ir 4t 

3(ai-a 2 ) 3(ai+a 2 ) 
4tt 8tt 



(22) 



where the aj = gf/(4:7i) and A is the Higgs self coupling defined by the potential V(4>) = 
A[(0V) - ^/ 2 ] 2 . 

Using this result for 7^ and the one-loop /3 functions for a±, (X2, and the lepton Yukawa 
couplings, we solve the RG equations to determine the operator coefficients C|(/-t) as a 
function of their values at the scale A. As in Ref. we find that the the running of the 
gauge and Yukawa couplings has a negligible impact on the evolution of the (/i). It is 
instructive to consider the results obtained by retaining only the leading logarithms ln(/x/A) 



5 The term in 733 proportional to A differs from that of Ref. 10], which contains an error. However, this 
change does not affect the bounds on the neutrino magnetic moments obtained in that work. 
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FIG. 3: One loop graphs for the mixing among n = 6 operators. Notation is as in previous figures. 
Various types of mixing (a-g) and self-renormalization (h-j) are as discussed in the text. 
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FIG. 4: Two-loop graphs for the mixing of the n = 6 operators. Only representee graphs for the 
mixing of the four-fermion operators 0^ ABCD into AD are shown. 



and terms at most first order in the Yukawa couplings. We find 
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F ABBD (A) 

3(«i - a 2 ) 

47T 



3(ai + a 2 ), /i 

1 In — 

8?r A 



F, BABD 



(A) In 



A 



+(9W2tt) [CS,AD(A)a 1 + CSr iU3 (A)a 2 ] ln| 



C" 



F, BABD 



(A) 



C 



_ 3(qi + a 2 ) ln /i 
8tt n A 

/ . ABBD (A) ln- 



3(ai - a 2 ) 

4:71 

■Wbb/Ax) [^ D (A)a 1 + ^ AD (A)a 2 ] In £ 



(23) 
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where 



C\ii) = a x C% A M - 3a 2 C 6 mA M (24) 

7± = (7i3 ±723) /2 
7 = 3(«i + 3«2)/167r 

We note that the combination of coefficients C+(v) enters the neutrino magnetic moment. 
Its RG evolution was obtained in Ref. [10] to zeroth order in the Yukawa couplings; here we 
obtain the corrections that are linear in f AA and Jbb- The corresponding contributions to 
the neutrino mass matrix 5m^ D and magnetic moment matrix jj,^ D are then given by 

AD _ ( 1,3 \ ri6 



S< U = - ) C b M>AD (v) (25) 

W2(^)Re{C»} . (26) 



a n 



From Eqs. ()23p25p26j) we observe that to linear order in the lepton Yukawa couplings, 

ad{^) receives contributions from the two magnetic moment operators and but not 
from the four fermion operators. This result is consistent with the result obtained by the 
authors of Ref. Q, who computed one-loop graphs containing the four- fermion operators 
of Eq. (JTJ) using massive charged leptons and found that contributions to m u oc m|. In 
the effective theory used here, the latter result corresponds to a one-loop computation with 
three insertions of the Yukawa interaction. However, mixing with was not considered in 

Ref. Q, and our result that this operator mixes with AD to linear order in the Yukawa 
couplings represents an important difference with the former analysis. 

We agree with the observation of Ref. [j| that the four fermion operators can mix with 

Ofj AD to linear order in the f AA via two- loop graphs, such as those indicated in Fig. HI 
These graphs were estimated in Ref. [y] by considering loops with massive W ± and Z° 
bosons that correspond in our framework to the diagrams of Fig. We observe, however, 
that the two-loop constraints will be weaker than those obtained by one-loop mixing with 
ad by ~ (a/ Aix)(v / K) 2 (modulo logarithmic corrections), so we do not consider this 
two-loop mixing in detail here. Moreover, because we work at a scale fi > v for which the 
use of massless fields is appropriate, and because we adopt a basis in which the Yukawa 
matrix and gauge interactions are flavor diagonal (but m^ D is not), the operators Op\ 12D 

and Op\ 21D will not mix with ad even a ^ two- loop order. 



IV. NEUTRINO MASS CONSTRAINTS 

To arrive at neutrino mass naturalness constraints on the g] coefficients, it is useful to 
tabulate their relationships with the dimension six operator coefficients. In some cases, one 
must perform a Fierz transformation in order to obtain the operator structures in Eq. (pj. 
Letting 

9l = -*Q 2 Ct (27) 
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we give in Table HP the ks corresponding to the various dimension six operators. 

Using the entries in Table HP and the estimates in Eqs. (fTtjj) . we illustrate how the bounds 
in Table H] were obtained. For the operator of l22D , for example, we have from Eqs. (j5H6|) 

l^,i 22 DUl6vr 2 (^) (28) 



5m}P \ f v\ 2 , T , 9 / 5m}, D \ ( v 



leading to 

where <5m AD denotes the radiative contribution to m^ D . Choosing A = v and 5m}, D = leV 
(corresponding to the scale of upper bounds derived from 3 H /3-decay studies ["ill. [15|) leads 
to the bounds in the first row of Table Similar arguments yield the other entries in the 
table. Note that the bounds become smaller as A is increased from v. 

The constraints on the g\ R RL that follow from mixing among the n = 6 operators follows 
straightforwardly from Eqs. (}2*3l l2*3j) and Table |Hj We obtain 



A similar expression holds for g\ L but with m M — > m e and 5m 2 u D — > 5w^P . Note that in 
arriving at Eq. (J3Uj) we have ignored the running of the AD {^) between A and v, since 

the impact on the Q^rrx, ls higher order in the gauge and Yukawa couplings. To derive 
numerical bounds on the gYnRi from Eq. ()30|) we use the running couplings in the MS 
scheme a = &(Mz) ~ 1/127.9, sin 2 9]y(Mz) ~ 0.2312 and the tree-level relation between 
the Higgs quartic coupling A, the Higgs mass m H , and v: 2A = (m H /v) 2 . We quote two 
results, corresponding to the direct search lower bound on itlh ~ 114 GeV and the one-sided 
95 % CL. upper bound from analysis of precision electroweak measurements, m# < 186 
GeV[u|. We obtain 

| 9 v I = fln-V' I 12 * m " = 1MGeV (31) 

|yLi?l V le V/V vj [7.5 xlO- 6 , m H = 186GeV V ; 

I y | f 5ml D \ / AV 1 |2.5 x 10" 4 , m H = 114 GeV 
\9rl\ 1 y leV J \ n v) |l.5xl0~ 3 , m H = 186GeV 

For A ~ 1 TeV, the logarithms are 0(1) so that for 5m v ~ 1 eV, the bounds on the Qlrrl 
derived from n = 6 mixing are comparable in magnitude to those estimated from mixing 
with the n = 4 mass operators. 

Although the four fermion operators do not mix with AD at linear order in the Yukawa 

couplings, they do contribute to the magnetic moment operators Og AD and 0$ AD at this 
order. From Eqs. (j23l I26|) we have 



~ g^2 (~) (x) R<3 [fAA C F,AAAD + }bB C F,ABBd\ ln ~ > ( 32 ) 
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TABLE II: Coefficients k that relate gl^ to the dimension six operator coefficients C\ via Eq. (|27|) , 



K 9lr 9lr 9rl 9rl 9lr 9rl 



C 6 F,X22D 1/4 1/8 - - - 

Cf,212D - - - - 

C 6 F ,n2D 3/4 1/8 - 

C F>211D - - 1/4 1/8 - 

^F, 121D " " 1/2 " 

C F ,221D - - 3/4 1/8 - 

C hn - - - - ~ 1 ' 2 - 

( v,n - 1 2 



where 5/j,^ d denotes the contribution to the magnetic moment matrix and /ig is a Bohr 
magneton. While O f \ aad does not contribute to /i-decay, the operator O f \ bbd does, and 
its presence in Eq. ()32|) implies constraints on its coefficient from current bounds on neutrino 
magnetic moments. The most stringent constraints arise for A — 1, B — 2 for which we find 

Current experimental bounds on \uf? p / Ur\ range from ~ 1CT 10 from observations of solar 
and reactor neutrinos 17, 0, 0, to ~ 3 x 10~ 12 from the non-observation of plasmon 
decay into vv in astrophysical objects j2l|. Assuming that the logarithm in Eq. (JSHj) is of 
order unity, these limits translate into bounds on g^ R and g\ R ranging from ~ 1 — * 0.03 
and ~ 0.3 — > 0.01, respectively. The solar and reactor neutrino limits on |/U^ xp //Ub| imply 
bounds on the g£ R that are weaker than those obtained from the global analysis of /i-decay 
measurements, while those associated with the astrophysical magnetic moment limits are 
comparable to the global values. Nevertheless, the bounds derived from neutrino magnetic 
moments are several orders of magnitude weaker than those derived from the scale of neutrino 
mass. 

The naturalness bounds on the C\ associated with the scale of m v have implications 
for the interpretation of /x-decay experiments. Because the coefficients C F U2D and C F 221D 

ST ' ' 

that contribute to g£ R R p are not directly constrained by m„, none of the eleven Michel 
parameters is directly constrained by neutrino mass alone. Instead, it is more relevant 
to compare the results of global analyses from which limits on the g 1 are obtained with 
the m u naturalness bounds, since the latter imply tiny values for the couplings gLR t RL- 
Should future experiments yield a value for either of these couplings that is considerably 
larger than our bounds in Table HI the new physics above A would have to exhibit either 
fine-tuning or a symmetry in order to evade unacceptably large contributions to m v . In 

S T 

addition, should future global analyses find evidence for non-zero g LRRL with magnitudes 
considerably larger than given by the m^-constrained contributions listed in Table U then 
one would have evidence for a non-trivial flavor structure in the new physics that allows 
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considerably larger effects from the operators O f ll2D and O f 221D than from the other four 
fermion operators. 

Finally, we note that one may use a combination of neutrino mass and direct studies of 
the Michel spectrum to derive bounds on a subset of the Michel parameters that are more 
stringent than one obtains from /^-decay experiments alone. To illustrate, we consider the 
parameters 5 and a, for which one has 

\ ~ P = \ \gU + \ \9ln\ 2 + §R* (9 S lh9 T lr) + (L~R) (34) 
a = 8 Re {g v RL (g s LR + Qg T L * R ) + (L <- R)} . (35) 

From Table HI we observe that the magnitudes of the g\ R RL contributions to p and a are 
constrained to be several orders of magnitude below the current experimental sensitivities, 

S T (6) (6) 

whereas the contributions g LRRL that arise from FU2D and F221D are only directly 
constrained by experiment. Thus, we may use the current experimental results for p to 
constrain the operator coefficients C F ll2D and C F 221D and subsequently employ the results 
- together with the m v bounds on the g\ RRL - to derive expectations for the magnitude 
of a. For simplicity, we consider only the contributions from C F U2D to p, and using the 
current experimental uncertainty in this parameter, we find 

In the parameter a, this coefficient interferes with 1£) : 

a = -6 (£f Re (c^C^ + ■ ■ ■) , (37) 

where the "+ • • • " indicate contributions from the other coefficients that we will assume to 
be zero for purposes of this discussion. From Eq. ()36|) and the m v limits on 1D we obtain 

H^x^GO 2 ®) ■ (38) 

For A = v, this expectation for |a| is more than two orders of magnitude below the present 
experimental sensitivity and will fall rapidly as A increases from v. A similar line of reasoning 
can be used to constrain the parameter a' in terms of m v and the CP-violating phases that 
may enter the effective operator coefficients. 



V. CONCLUSIONS 

The existence of the small, non-zero masses of neutrinos have provided our first direct 
evidence for physics beyond the minimal Standard Model, and the incorporation of m v 
into SM extensions is a key element of beyond-the-SM model building. At the same time, 
the existence of non-vanishing neutrino mass - together with its scale - have important 
consequences for the properties of neutrinos and their interactions that can be delineated 
in a model-independent manner 

0, M, H 0- In this paper, we have analyzed those 



implications for the decay of muons, using the effective field theory approach of Ref. [1 
and concentrating on the case of Dirac neutrinos. We have derived model-independent 
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naturalness bounds on the contributions to the Michel parameters from various n = 6 
operators that also contribute to the neutrino mass matrix via radiative corrections. 

Our work has been motivated by the ideas in Ref. 0, but our conclusions differ in im- 
portant respects. Importantly, we find — after properly taking into account SU(2) L xU(l)y 
gauge invariance and mixing between n = 6 and n = 4 operators that cannot be stud- 
ied using dimensional regularization — the dominant constraints on the contributions to 
the Michel parameters occur at one-loop order, rather than through two-loop effects as in 
Ref. 0. Consequently, the bounds we derive are generally two orders of magnitude (or 
more) stronger than those of Ref. |9j]. In addition, we carefully study the flavor structure 
of the operators that can contribute to /i-decay and and find that there exist four-fermion 
operators that do not contribute to the neutrino mass matrix through radiative corrections. 

S T /ITT 

These operators contribute to the effective scalar and tensor couplings g LR RL of Eq. . In 
contrast, all operators that generate the g\n rl terms contribute to m^ D , so these effective 
couplings do have neutrino-mass naturalness bounds. From a model-building perspective 
it might seem reasonable to expect the coefficients of the unconstrained four-fermion oper- 
ator coefficients to have the same magnitude as those that are constrained by m u , but is 
important for precise muon-decay experiments to test this expectation. 

While we have focused on the implications of Dirac mass terms, a similar analysis for 
the Majorana neutrinos is clearly called for. Indeed, in the case of neutrino magnetic mo- 
ments, the requirement of flavor non-diagonality for Majorana magnetic moments leads to 
substantially weaker naturalness bounds than for Dirac moments. While we do not antici- 
pate similar differences between the Majorana and Dirac case for operators that contribute 
to /i-decay, a detailed comparison will appear in a forthcoming publication. 
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